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ABSTRACT 

Cassiopeia A the youngest supernova remnant known in the Milky Way is one of the brightest 
radio sources in the sky and a unique laboratory for supernova physics. Although its compact 
remnant was discovered in 1999 by the Chandra X-Ray Observatory, nowadays it is widely 
accepted that a neutron star lies in the center of this supernova remnant. In addition, new ob¬ 
servations suggest that such neutron star with a low magnetic field and evidence of a carbon 
atmosphere could have suffered a hypercritical accretion phase seconds after the explosion. 
Considering this hypercritical accretion episode, we compute the neutrino cooling effect, the 
number of events and neutrino flavor ratios expected on Hyper-Kamiokande Experiment. The 
neutrino cooling effect (the emissivity and luminosity of neutrinos) is obtained through nu¬ 
merical simulations performed in a customized version of the FLASH code. Based on these 
simulations, we forecast that the number of events expected on the Hyper-Kamiokande Ex¬ 
periment is around 3195. Similarly, we estimate the neutrino flavor ratios to be detected con¬ 
sidering the neutrino effective potential due to the thermal and magnetized plasma and thanks 
to the envelope of the star. It is worth noting that our estimates correspond to the only trust¬ 
worthy method for verifying the hypercritical phase and although this episode took place 330 
years ago, at present supernova remnants with these similarities might occur thus confirming 
our predictions for this phase. 
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1 INTRODUCTION 

At a distance of d z — 3.4lo'i kpc, Cassiopeia A (Cas A) is one 
of the closest and youngest known supernova remnants (Reed et al. 
1995; Chang & Bildsten 2004). It was related to the supernova ob¬ 
served by John Flamsteed in 1680 (Ashworth 1980). As was re¬ 
ported by Krause et al. (2008), Cas A supernova remnant was pro¬ 
duced in a lib type supernova explosion. The progenitors of such 
explosions, as was showed by Chevalier (2005), are red supergiants 
which already have lost their hydrogen envelope via a powerful 
stellar wind. Then, it is inferred that such progenitors with masses 
of the order of 20 M 0 must have a helium core of 6 M 0 and an iron 
core of 1.45 M 0 . This scenario should produce a neutron star with 
1.4 M 0 . It seems to be the case of Cas A. 

Studying the atmosphere surrounding the neutron star of Cas A, Ho 
& Heinke (2009) found evidence of carbon in the composition of 
this atmosphere. It suggests that a strong accretion could have oc¬ 
curred after the neutron star formation. In the core-collapse frame¬ 
work other evidences related to the low strength of magnetic field, 
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i) the non-detection of the spectral features related to the electron 
cyclotron resonance in the Chandra energy range (which is usu¬ 
ally associated to a B ~ (1 — 5) x 10 11 G), ii) the absence of a 
visible pulsar wind nebula and iii) no samples of magnetospheric 
activities (such as radio or 7-ray fluxes) from the direction of Cas 
A, as found in pulsars with B ^ 10 12 G, suggest that the strong 
magnetic field was buried because of a strong accretion (Muslimov 
& Page 1995). If the core-collapse paradigm is correct, the shock 
wave sweeps the outer layers of the progenitor until it encounters a 
discontinuity in density. At this point, a reverse shock is generated 
leading to a fallback episode which allows to deposit large amounts 
of material on the newborn neutron star surface. As a consequence 
of the fallback episode, the magnetic field is submerged within the 
new crust formed during such episode (Bernal et al. 2013). 

In the hypercritical phase, super-Eddington accretion gives rise to 
the material of core-collapse to become dense enough so that pho¬ 
tons are trapped. Therefore, the gravitational energy is carried out 
due to neutrino emission in the formation of the newborn neu¬ 
tron star. When the accretion rate M exceeds by several orders 
of magnitude the Eddington accretion rate Msdd 10 18 gr s _1 
so-called ’’Eddington limit”, then it is denoted as hypercritical ac¬ 
cretion (Blondin 1986), and the accretion rate required is higher 
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than ~ 10 3 M E dd■ 

In the hypercritical accretion episode (seconds after the supernova 
explosion), due to the inverse beta decay, electron-positron an¬ 
nihilation ( e~ + e + Z i/j + Pj) and nucleon-nucleon 
bremsstrahlung (N + N —>N + N + i/j + 9j)forj. = e, i/, r, ther¬ 
mal neutrinos (with energies in the range 5 MeV^ E v ^ 30 MeV) 
will be produced at the core and then oscillate in their propagation 
through the collapsing star. The resonant conversion (from one fla¬ 
vor to another) and the properties of these neutrinos will get modi¬ 
fied due to neutrino effective potential generated by the magnetized 
and thermal plasma (Sahu et al. 2009b, a), and later by the collaps¬ 
ing and surrounding material of the progenitor. It is worth noting 
that, although neutrino cannot couple directly to the magnetic field, 
its effect can be felt through the coupling to charge particles in the 
thermal plasma. Therefore, the values of neutrino flavor ratios cre¬ 
ated at the core will be different from the values on the surface of 
the star and on Earth. 

Recently, Fraija et al. (2014) revisited Chevalier’s model in a nu¬ 
merical framework to focus on the neutrino cooling effect on the 
SN1987A fail-back dynamics. They estimated the number of neu¬ 
trinos that must have been seen on the Super-Kamiokande neutrino 
experiment, hours after the ~ 20 events detected from the direc¬ 
tion of SN1987A. In the current work, we consider that Cas A 
supernova remnant underwent the hypercritical accretion episode 
in its early life as neutron star. We present a numerical study of 
hypercritical accretion of matter onto the neutron star surface. As 
an evidence of this episode we estimate the neutrino luminosity, 
the expected number of neutrinos and the neutrinos flavor ratio on 
Earth. We employ the FLASH code to carry out the numerical sim¬ 
ulations of the reverse shock and the complex dynamics near the 
stellar surface, including neutrino cooling processes, a more de¬ 
tailed equation of state, and an additional degree of freedom in the 
system. The paper is arranged as follows. In section 2 we describe 
the dynamics of the hypercritical accretion phase. In section 3 we 
study the thermal neutrino expectation. In section 4 we discuss the 
results and a brief conclusions is drawn in section 5. 


2.1 The new crust on the neutron star surface 


When the hypercritical phase is occurring, the piled material 
builds a small region at the same height scale where the neutrino 
cooling takes place. It forms a new thin crust with a strong magnetic 
field immersed within it. In this crust formed at a region r ns ^ 
r ^ ( r ns + r c ) on the surface of the newborn neutron star surface, 
photon field radiation and neutrinos are thermalized to a value in 
the range 1 MeV ^ T ^ 5 MeV. Here, r ns = 10 6 cm is the 
neutron star radius and r c ~ 200 — 400 m is the size of the new 
thin crust. The energy density and neutrino emissivity are given by 
(Dicus 1972) 


U b (T) = 3.77 x 10 26 
and 

i v = 0.97 x 10 2 


MeV 


erg cm 


MeV 


-1 -3 

erg s cm , 


( 1 ) 

( 2 ) 


respectively. At a first approach and due to the high dependence of 
the pair annihilation process on pressure near the stellar surface, 
the neutrino emissivity can be written as 

e„ = 1.83 x 10 p erg cm s , (3) 


where the total pressure p includes (besides the gas pressure) the 
pressure due to electron/positron pairs which is defined by 



(4) 


with a the radiation constant. It is worth noting that, in the hyper¬ 
critical phase, the magnetic field may be submerged and confined 
into the new crust (formed by the piled material). The strength of 
the magnetic field could reach values in the range of 10 11 G ^ 
B < 10 13 G. Also, the density and pressure of radial profiles can 
be described as power laws, whereas the mean velocity is null in 
this region. 


2 DYNAMICS OF THE HYPERCRITICAL ACCRETION 
PHASE 

In the framework of the stellar core-collapse, a reverse shock is 
generated when the expanding shock wave encounters dense outer 
layers on its path. Whereas the forward shock advances out of the 
star, the reverse shock leads to a fall-back phase which can induce 
the hypercritical accretion onto the newborn neutron star surface 
during the following 3 hours after the explosion. This scenario is 
only possible if the progenitor has a tenuous H/He envelope sur¬ 
rounding a dense He core, as suggested by Ho & Heinke (2009). In 
this approach, the neutrino emission determines the formation of a 
quasi-hydrostatic envelope around the newborn neutron star in the 
hypercritical phase. In the above mentioned scenario, we can iden¬ 
tify four dynamically different regions: (i) a thin cmst surrounding 
the newborn neutron star, which is formed during the hypercritical 
phase; (ii) an envelope in quasi-hydrostatic equilibrium formed by 
the third expansive shock; (iii) a free-fall region above this enve¬ 
lope; and (iv) the external layers of the progenitor. Now, we are 
going to consider the dynamics inside each region around the com¬ 
pact remnant. 


2.2 Quasi-hydrostatic envelope 


The reverse shock induces the hypercritical accretion onto the 
newborn neutron star surface through the fall-back episode. This 
accreted material will bounce against the surface of the newly 
born neutron star, building a new expansive shock. This expansive 
shock builds an envelope in quasi-hydro static equilibrium, with 
free falling material raining over it. The structure of such envelope 
in quasi-hydro static equilibrium is given by (Chevalier 1989) 



where the subscript 5 refers to the value of density (p s ), pressure 
(p s ) and velocity (v s ) at the shock front. Given the typical values 
of the neutron star: a mass M — 1.4 M 0 , an accretion rate M — 
500 Mq yr -1 and a radius r ns = 10 6 cm, the shock radius is 


r s 


7.73 x 10 8 cm 


M 

1.4 Mq 


-0.04 


/ r ns \ 1-48 
V10 6 cm/ 


X 


M 


500 M 0 yr - 


-0.37 


•( 6 ) 


The pressure p s and density p s are determined by the strong shock 
condition p s — f po^o and ps = 7po, respectively and v s by the 
mass conservation v s = — \vq. Here po and vo are the density and 
velocity just outside the shock, respectively. In this case, the radial 
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distance of the accretion shock that is controlled by the energy bal¬ 
ance between the accretion power and the integrated neutrino losses 
is given by 

GMM w 

— R — = J ev ^ r h lr - C) 

Chevalier (1989) approached the size of the envelope as R /4 for the 
pressure profile in an atmosphere in quasi-hydro static equilibrium. 
With this, and using the analytical formula of Dicus (1972) for the 
neutrino emissivity, the shock radius can be calculated through the 
energy balance as 


GMM 

R 


= 4t tR 2 



( 8 ) 


Then, the high pressure near the neutron star surface becomes 
Pns — 1-86 x 10 _12 dyncm _2 Mrf //2 . (9) 


It allows the pair neutrino process to be the dominant mechanism 
in the neutrino cooling. From the strong shock conditions and eq. 
(5), we get 

pqhe(r) = 7.7 x 10 2 ^) gem -3 , (10) 

where the quasi-hydrostatic envelope radius lies in the range ( r ns + 

r c ) < r < r s . 


2.3 Free fall zone 

In this case, material at r s ^ r ^ 77 begins to fall with the velocity 
and density profiles given by 

= and Pff (r) = —M^. (11) 

Taking the typical values for such object, M — 1.4 M© and M — 
500 M© yr -1 , the velocity is written as 

v ff( r ) ~ 7.81 x 10 7 ^ — ^ ems -1 , (12) 

and the density of material in free fall as 

Pff{r) = 5.74 xlO -2 ^^ gem -3 , (13) 

where 77 — 6.3 x 10 10 cm. 


2.4 External layers 


By considering a supergiant profile (typical for a pre-supernova), 
the density of the outer layers can be described as (Chevalier & 
Soker 1989) 


Pelf) 


3.4 x 10 5 gcm 3 

x f(-r; 

x \ ( R*\ n 1 (r-i?*) n 2 _ 
( l r ) (r b -R *) n 2 > 


rh < r < n 
r > rb. 


’(14) 


with rb — 10 12 cm, i?* ~3 x 10 12 cm, m — 17/7 and n 2 — 5. 


3 THERMAL NEUTRINO EXPECTATION 

During the hypercritical accretion episode a copious amount of 
neutrinos in the energy range of 1 - 30 MeV is produced. The pro¬ 
cesses of neutrino emission involved in the episode are: 

• Electron-positron annihilation process, pairs are annihi¬ 
lated to create neutrino-antineutrino pairs ( e~ + e + — > v + z7). 

• Photo-neutrino process. A photon is Compton scattered to pro¬ 
duce a neutrino-antineutrino pair (7 + —>> e ± + z/ + z/). 

• Plasmon decay process. A photon propagating within an 
electron plasma (plasmon) is spontaneously transformed into a 
neutrino-antineutrino pair (7 —>> v + v). 

• Bremsstrahlung process. A neutrino-antineutrino pair is cre¬ 
ated by either an electron-nucleon interaction (e^ + N —y e ± +A^+ 
v + v) or nucleon-nucleon interaction (N + N ^ N + N + V + V). 


3.1 Number of Expected Neutrinos 


The expected number of neutrinos can be estimated requir¬ 
ing the neutrino luminosity of the hypercritical phase L — 
4tv d 2 z E 2 dN/dE and the effective volume of the neutrino detector 
V. The expected event rate can be written as 


N ev = VN A p N jj alf^dEdT, (15) 

where Na = 6.022 x 10 23 g -1 is the Avogadro’s number, pN = 
2/18gem -3 is the nucleons density in water (Mohapatra & Pal 
2004), crcc P — 9 x 10 -44 El e /MeV 2 is the cross section (Bahcall 
1989) and dT is the detection time of neutrinos. The number of 
events (eq. 15) is 


N ev 


< > 
t 


47 \d 2 z < E 9p > 


VNAPNCTcf <E V 
VNaPncjI? 


> 2 < 


dN 

dE Vp 


(16) 


where we have assumed an average energy < E„ e E Ve during 
a period of time t (Mohapatra & Pal 2004). 


3.2 The Hyper-Kamiokande Experiment 

The Hyper-Kamiokande detector will be the third generation wa¬ 
ter Cherenkov in Kamioka, designed for a vast variety of neutrino 
studies. At 8 km south of its predecessor Super-Kamiokande, the 
Hyper-Kamiokande detector will be located in the Tochibora mine 
of the Kamioka Mining and Smelting Company, near Kamioka 
town in the Gifu Prefecture, Japan. This detector consists of two 
separate caverns, each having an egg-shaped cross section of 48 
meters wide, 54 meters tall, and 250 meters long. The entire array 
consists of 99.000 PMTs, uniformly surrounding the region. It will 
have a total (fiducial) mass of 0.99 (0.56) million metric tons, ap¬ 
proximately 20 (25) times larger than that of Super-Kamiokande. 
Among the physical potentials of this detector is the detection of 
astrophysical neutrinos from the Supernova remnant and the stud¬ 
ies of neutrino oscillation parameters (Abe & et al. 2011a; Hyper- 
Kamiokande Working Group et al. 2014). 


3.3 Neutrino Oscillation 

The properties of these neutrinos are modified when they propa¬ 
gate in matter. In the following subsections we will describe the 
dynamics of neutrinos going through matter. 
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3.3.1 Two-Neutrino Mixing 


The neutrino oscillations for two-neutrino mixing ( v e «<-)> z/ MjT ) are 
determined by the evolution equation in matter given by 

*( v *) =-i 5 IH-( 4 N^ L - 2 cos2e sin20\^ey (1?) 
AE U \ sin 20 0 J J 


where 14 // is the neutrino effective potential generated by the 
medium, E v is the neutrino energy, 5m 2 is the neutrino mass dif¬ 
ference and 6 is the neutrino mixing angle. From the conversion 
probability P Vc ^ v ^ Vr ){t) = sin(26») sin (wt/ 2)) 2 with 

u = If- \J ( 2 4n^ ~ ~ cos 2 #) 2 + si 1 ' 2 29, the oscillation length 
of neutrino is written as 

Lose = . = =, (18) 

V ( 2E 6mi f ' L - c °s20 ) 2 + sin 2 29 

where L v = 4-KEo/Sm 2 is the oscillation length in vacuum. The 
resonance condition and resonance length can be written as 


Veff 

and 

Lres 


5m 2 

2K 


cos 2 6 , 


L v 

sin 20 ’ 


(19) 


( 20 ) 


respectively. 


3.3.2 Three-neutrino Mixing 


The dynamics of neutrino oscillations in matter for three-neutrino 
mixing is determined by 

^ = -i[U ■ Hq ■ U f + diag(V ef f, 0,0 )]v, ( 21 ) 

where Hq = ^-diag{—8m |i, 0 , 5m 2 2 ), U is the matrix of three- 
neutrino mixing (Gonzalez-Garcia & Nir 2003; Akhmedov et al. 
2004; Gonzalez-Garcia & Maltoni 2008; Gonzalez-Garcia 2011) 
and the state vector in the flavor basis is V = (v e ,v^,is r ) T . The 
survival and conversion probabilities for electron, muon and tau 
are given in Fraija (2014). The oscillation length of neutrino for 
three-neutrino mixing is given by 


Lose - 


^/cos 2 2 fl 13 (l ~ ) 2 + sin2 2 ^ 13 ’ 


( 22 ) 


where L v = 4:7rE u /8m 2 2 . The resonance condition and resonance 
length are given by 


i- F 2 T vVef L = o, 


and 

L r es — 


8m^ 2 cos 2^13 
L v 


sin 2^13 ’ 


(23) 


(24) 


respectively. 


3.3.3 Neutrino Oscillation from Source to Earth 

Neutrinos oscillate in their paths to Earth. In this case, the dy¬ 
namics is governed by the Schrodinger equation in vacuum and 
the Hamiltonian is given as a function of neutrino mixing param¬ 
eters. The conversion probability from a flavor state a to a fla¬ 
vor state (3 is written as = |< i/p(t)\i/ a (t = 0 ) >| = 


S a p - UaiUpiUajUpi sin 2 ( SmLL/AEr^j, where Uji is 

once again the three neutrino mixing (Gonzalez-Garcia & Nir 
2003; Akhmedov et al. 2004; Gonzalez-Garcia & Maltoni 2008; 
Gonzalez-Garcia 2011). Here the sin term in the probability has 
been averaged to ~ 0.5 for distances (L) longer than the Solar Sys¬ 
tem (Learned & Pakvasa 1995). 


4 RESULTS 


4.1 Plasma-dynamics on the stellar surface 


We carry out hydrodynamic numerical simulations of the hypercrit¬ 
ical accretion phase in a two-dimensional spherical mesh, using the 
Flash method developed by Fryxell et al. (2000). This multi-physics 
Eulerian and parallelized code was designed to figure out many 
problems in high-energy astrophysics, especially those related to 
the dynamics of plasmas around neutron stars. Recently, using a 
customized version of such numerical code, Fraija et al. (2014) pre¬ 
sented numerical simulations of strong matter accretion onto the 
newborn neutron star to compute the size of the neutrino-sphere, 
the emissivity and the luminosity of neutrinos in the SN1987A sce¬ 
nario. Requiring adequate parameters to describe the Cas A sce¬ 
nario, we will follow a similar procedure to that performed for 
SN1987A. 

The Flash code uses the piecewise-parabolic method (PPM) solver 
to solve the whole set of hydrodynamic equations for a com¬ 
pressible gas dynamics in one, two, or three spatial dimensions. 
This directionally-split method makes use of a second-order Strang 
time splitting, which uses a finite-volume spatial discretization of 
the Euler equations together with an explicit forward time differ¬ 
ence. Time-advanced fluxes at cell boundaries are computed using 
the numerical solution to a Riemann shock-tube problem in each 
boundary. Initial conditions for each Riemann problem are deter¬ 
mined by assuming the non-advanced solution to be piecewise- 
constant in each cell. Requiring the Riemann solution, the effect 
of explicit nonlinearity is introduced into the differential equations. 
It allows the calculation of sharp shock fronts and contact disconti¬ 
nuities without introducing significant nonphysical oscillations into 
the flow. Since the value of each variable in each cell is assumed 
to be constant, this method is limited to first-order accuracy in both 
space and time. PPM improves on Godunovs method by represent¬ 
ing the flow variables with piecewise-parabolic functions. It also 
uses a monotonicity constraint rather than artificial viscosity to con¬ 
trol the oscillations near discontinuities. The Euler equations can be 
written in conservative form as 


dpw 

dt 

dpE 

dt 


ft +V - ( '’V) 

= 0 

(25) 

+ V • (pvv) + VP 

= pg 

(26) 

+ V • [v (pE + P)] 

= pv • g , 

(27) 


where the previous equations (25, 26 and 27) represent the conser¬ 
vation of mass, momentum and total energy, respectively. In this 
context, p is the fluid density, v is the fluid velocity, P is the pres¬ 
sure, g is the acceleration due to gravity, t is the time coordinate 
and E represents the sum of the internal energy e and kinetic en¬ 
ergy per unit mass 


E = z+\ |v| 2 . 


(28) 


The pressure is obtained from the energy and density using the 
equation of state. In this case, we use a customized version of 


© 0000 RAS, MNRAS 000, 000-000 















Hypercritical accretion phase and neutrino expectation in the evolution of Cassiopeia A 5 


the Helmholtz Unit of the Flash code which includes contributions 
from the nuclei, pairs, radiation and the Coulomb correction. 
This routine is appropriate for addressing astrophysical phenomena 
in which electrons and positrons may be relativistic and/or degen¬ 
erated, and the radiation may significantly contribute to the ther¬ 
modynamic state. The pressure and internal energy are calculated 
considering all components 

Ptot — Prad + Pion H - Pele H - Ppos + Pcoui , (29) 

ttot — ^rad + tion tele H - tpos H - t cou l . (30) 


Here the subscripts “rad”, “ion”, “ele”, “pos” and “coul” represent 
the contribution from radiation, nuclei, electrons, positrons, and 
corrections for Coulomb effects, respectively. The radiation por¬ 
tion assumes a blackbody in local thermodynamic equilibrium, the 
ion portion (nuclei) is treated as an ideal gas with 7 = 5/3, and the 
electrons and positrons are treated as a noninteracting Fermi gas. 
That is, the number densities of free electrons N e i e and positrons 
N pos in the noninteracting Fermi gas formalism are given by 

/2 

N e le = —2 m 3 el3 3/2 [F 1 / 2 (Vele, P) + F 3 / 2 (r)ele, /?)] l (31) 

Npos = ^rn 3 e p 3/2 [Ji/afapo.,0) + PF 3/2 ( Vpos , P)] , (32) 

where m e is the electron rest mass, ft = T/m e is the normalized 
temperature, rj e i e = p/T is the normalized chemical potential for 
electrons, ?7p OS = —r] e i e —2//3 is the normalized chemical potential 
for positrons, and Fk(p, /3) is the Fermi-Dirac integral given by 

00 

r ■ 

0 

Because the electron rest mass is not included in the chemical po¬ 
tential, the positron chemical potential must have the form given by 
if pos • For complete ionization, the number density of free electrons 
in the matter is 


N e le, matter — — -N a p — Z Ni 0 n , 

A 


and charge neutrality requires 


(34) 


N et 


e,matter 


— N e l e Np OS • 


(35) 


Solving the previous equations with a standard one-dimensional 
root-finding algorithm we determine 77, and the Fermi-Dirac inte¬ 
grals can be evaluated given the pressure, specific thermal energy, 
and entropy due to the free electrons and positrons. Full details of 
the Helmholtz equation of state are provided in Timmes & Swesty 
( 2000 ). 

As was pointed out, neutrino energy losses are dominated by the 
annihilation processes which involve the formation of neutrino- 
antineutrino pairs through the annihilation near the stellar sur¬ 
face (e + + e~ —>• v + v). In addition, we include other relevant 
neutrino processes present in such regime, which are implemented 
in a customized module and described in Itoh et al. (1996). We ig¬ 
nore neutrino absorption and heating, and nuclear reactions. 
Recently, Bernal et al. (2010, 2013) showed that, regardless of its 
initial configuration, the magnetic field is buried and submerged 
under the stellar surface during the hypercritical phase. Therefore, 
as the magnetic field is submerged within the surface, it does not 
play an important role in the dynamics of quasi-hydro static enve¬ 
lope, although the dynamics itself of the envelope is relevant. Here, 
we do not take into account the magnetic field and just focus on the 


hydrodynamics near the stellar surface. The hydrodynamic numeri¬ 
cal simulations of the hypercritical accretion regime are carried out 
in a two-dimensional spherical mesh (r, 0 ), where the radial com¬ 
ponent r lies in the range 10 6 cm ^ r ^ 3 x 10 6 cm. Due to 
the symmetry of the problem, we only simulate a quarter of the to¬ 
tal domain (7 t/4 ^ 0 ^ 37 t/4 ) which has 2048 x 2048 effective 
zones as spatial resolution. It is worth noting that it corresponds 
to the maximum level resolution which is equal to 6 for this work. 
The time resolution of these simulations is dt ~ 10 -7 s. We use, 
for the discretization grids, a block-stmctured oct-tree based on an 
adaptive grid used by the PARAMESH library which is included 
in FLASH code. If the Adaptive Mesh Refinement (AMR) grid is 
used, then the formation of the physical domain starts at the lowest 
level of refinement. Initial conditions are applied to each block at 
this level by calling initial condition routines. The Grid Unit then 
checks the refinement criteria in the blocks that it has created and 
also refines these blocks when the criteria are met. It then calls 
the routines of the initial conditions to initialize the newly created 
blocks. This process repeats until the grid reaches the required re¬ 
finement level in the areas marked for refinement. 

In the boundary conditions, we impose mass inflow along the top 
edge of the computational domain. We assume that the accreting 
matter is not magnetized. Concerning the FLASH code, the parallel 
structure of blocks means that each processor works independently. 
If a block is on a physical boundary, the guard cells are filled by cal¬ 
culations since there are not any neighboring blocks from which to 
copy values. At the top of the domain, initially we set the velocity in 
free fall to be v r = 7/ 2GM/r in all the guard cells, and set a den¬ 
sity profile which fixes a constant accretion rate, p — M/Anr 2 v r , 
as in the computational domain. The temperature is uniform and 
set to T = 10 9 K. Other variables are calculated with the equa¬ 
tion of state. When the shock reaches this boundary, we allow it to 
leave the computational domain. At the same time, this boundary 
switches to the Quasi-hydrostatic Envelope mode, i.e, it is adapt¬ 
able for consistency with the model. Lor this purpose we use the p, 
p, and v profiles given by the analytical solution (eqs. 5 and 6). It 
is justified due to the strong initial transient and short duration of 
the hypercritical phase. At the bottom, on the neutron star surface, 
we use a custom boundary condition that enforces the hydrostatic 
equilibrium. In order to establish this boundary in the problem of 
hypercritical accretion, we fix the velocities as null in all the guard 
cells, (v r = vq = 0), and copy the density and the pressure of the 
first zone of the numerical domain. This zone corresponds to the 
neutron star surface (p = p(ic),p — p(ic) + pv 2 + pgh ), where 
ic is the first zone in the domain. The rest of thermodynamics vari¬ 
ables are calculated from the equation of state. In addition, we have 
implemented periodic conditions along the sides. Periodic (wrap¬ 
around) boundary conditions are initially configured in this routine 
as well. If periodic boundary conditions are set in the x-direction, 
for instance, the first blocks in the x-direction are set to have as their 
left-most neighbor the blocks that are the last in the x-direction, and 
vice versa. Thus, when the guard cell filling is performed, the peri¬ 
odic boundary conditions are automatically maintained. 

We focus our study on the dynamics of the initial transient when the 
reverse shock rebounds off the stellar surface and builds a quasi¬ 
hydrostatic equilibrium envelope, instead of following the shock. 
We start the simulation requiring the free-fall profiles described in 
subsection 2.3 and considering an initial temperature of 10 9 K. The 
initial value of temperature was elected because it is close to the 
value expected from the theoretical model, and then the code finds 
the correct pressure profile after some steps of simulation. The neu¬ 
tron star mass and the accretion rate, for the Cas A parameters, are 
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chosen as M ~ 1.4 M© and M — 500 M©yr _1 , respectively. 
Fig. 1 shows the time evolution color maps of density (A) and neu¬ 
trino emissivity (B) for different times: t = 0.2 ms (up) shows the 
initial transient of the third expansive shock with material falling 
onto it; t = 1 ms (right) illustrates the evolution of the shock in the 
computational domain resulting in a complex morphology above 
the stellar surface; t — 5 ms (down) exhibits the moment when the 
shock reaches the external boundary and leaves the computational 
domain; and finally, t = 30 ms (left) shows the quasi-hydrostatic 
envelope established with a new thin crust on the neutron star sur¬ 
face. The hydrodynamic instabilities and the rich morphology ob¬ 
served inside the envelope in early times have disappeared. Note 
that in the new thin crust the neutrino cooling is very effective cre¬ 
ating an energy sink that allows the material to be deposited on the 
surface slowly. That is, in such region the neutrino emissivity is 
highly efficient. 

In this analysis, not only hydrodynamic simulations of the hy¬ 
percritical regime are performed, but also the results of previous 
MHD simulations for similar parameters of Cas A are presented. 
Such results show a complete magnetic field submergence under 
the new formed crust, as it was reported in Bernal et al. (2010, 
2013). In these papers, to perform the MHD simulations the au¬ 
thors used the Split Eight-Wave solver to solve the whole set of 
MHD equations. They considered wide accretion columns in 2D 
and 3D to follow the accretion shock and also a magnetic field 
loop (in the shape of an hemi-torus) as magnetic initial condi¬ 
tion. On the central hemi-circle of the loop the field had a strength 
Bo = 10 12 G, it was shaped as a Gaussian, i.e., with strength 
B(d) = Bo exp [— (cI/Rl) 2 )], d being the distance to the loop 
central hemi-circle and Rl — 1 km. The two feet of the loop were 
centered at x — —5 and x — +5 km, and z = 0 in the 3D model. 
When magnetic field was present, it was put at the bottom bound¬ 
ary in such a form that it was continuous from the guard cell to the 
physical domain, i.e, the authors anchored the magnetic field onto 
the neutron star surface and in the rest of the guard cells it was null. 
On the other hand, Figure 2 shows radial profiles of density and 
temperature as a function of the stellar radius when the envelope 
has been well established. It is worth noting that density profiles ob¬ 
tained by MHD and HD are very similar, hence we show only those 
densities and temperatures obtained by MHD. This is justified be¬ 
cause in this case we were able to follow the evolution of the shock 
until it formed the quasi-hydrostatic envelope with matter falling 
over it. In the left-hand panels we show the density and temperature 
profiles obtained from simulations for several hyperaccretion rates, 
in terms of a fiducial accretion rate (M ~ 500 M 0 yr -1 ), whereas 
in the right-hand panels we show the analytical solutions (with the 
same rates), using the analytical model (Chevalier 1989). Note the 
excellent agreement between them. The three regions (new crust, 
envelope and free fall) are visible and distinct. In this case, label 
[1] represents the Cas A accretion rate, and labels [10] and [100] 
are accretion rates one and two orders of magnitude larger, respec¬ 
tively. The envelope is well established in 600 ms for label [1], 300 
ms for label [10] and 100 ms for label [100]. A detail comparison 
between the MHD and HD solvers was reported in Bernal et al. 
(2010). In the new crust zone, the magnetic field is compressed and 
amplified, as shown in Figure 3 (upper-left panel). The magnetic 
field submergence into the new crust formed in the hypercritical 
regime allows us to justify the fact that the magnetic field does 
not play an important role in the subsequent dynamics of the enve¬ 
lope. Both in the hydrodynamic and MHD cases the estimated size 
of the neutrino-sphere (including all the relevant neutrino cooling 
processes) is r ~ 3.7 x 10 5 cm ~ (1/3)7?. The mean value of 


emissivity in this region is e v ~ 2.4 x 10 30 ergs 1 cm 3 . The 
volume of the neutrino-sphere is given by V ~ (4/3)ttR 3 = 
(4/3)7t X 10 18 cm 3 and then the neutrino luminosity is obtained 
to be L u ~ 8.4 x 10 48 ergs -1 . This estimate agrees with the 
numerical value obtained by integrating directly the neutrino lumi¬ 
nosity in the whole computational domain. The upper-right panel 
in Figure. 3 shows the parameter-space of temperature and density. 
Taking into account the Itoh et al. (1996) tables, this panel has been 
divided into grey zones, highlighting the different regions where 
neutrino cooling processes dominate. Also, we have drawn lines in 
white color to show the region of interest for the Cas A parame¬ 
ters. The panel below in Figure 3 shows the time evolution of the 
neutrino luminosity integrated in the whole computational domain 
for the hydrodynamic case. In the initial transient the luminosity 
reaches a maximum value, showing then small oscillations about a 
fixed value. Similar results are found in the MHD case. 


4.2 Neutrino expectation 

Taking into account the values of neutrino Luminosity Lp e — 
8.4 x 10 48 erg s -1 , effective volume V ~ 0.56xl0 12 cm 3 (Hyper- 
Kamiokande Working Group et al. 2014) and the average neutrino 
energy < Ep e >~ 15MeV, from eq. (16) we get that the num¬ 
ber of events expected from the hypercritical accretion episode on 
Hyper-Kamiokande is 3195. 

These thermal neutrinos will oscillate when they propagate through 
the star. On the neutron star surface, the plasma is highly magne¬ 
tized with B~3x 10 12 G(seefig.3 left-hand panel) and thermal- 
ized at T ~ 3 MeV (see fig. 2 panels below). The neutrino effective 
potential at the moderate field limit is given by (Fraija 2014) 


Veff = 


\/2 Gp m g B 
^BBc 


£(-!)'sinhaJ -f- l + 4-| )K 1 (a l ) 


E 2 


^ An ( 2 + ^t( 3_2 ^ + 4 S)) Xl(<T ' An) } 


— 4—y- —— 1)* cosh a; / ^-Ko(rji) + ^ \ 2 n K 0 (<Ji\ n ) 


where K* is the modified Bessel function of integral order i, n is 
the chemical potential, A 2 = 1 + 2nB/B c , ai = (l + 1) fi/T, 
al — (l + 1) rrie/T and B c = 4.4 x 10 13 G. As shown in Fig¬ 
ure 4, the effective potential is an increasing function of magnetic 
field for T=1 MeV, 2.5 MeV and 5 MeV. The value of the effec¬ 
tive potential lies in the ranges: 2.3 x 10 -11 eV - 5.1 x 10 -lo eV 
for T=1 MeV, 2.61 x 10 -11 eV - 5.7 x 10 -lo eV for 2.5 TeV, and 
2.64 x 10 -11 eV - 6.3 x 10 -lo eV for 5 MeV when the magnetic 
field lies at 4 x 10 9 G ^ B ^ 4 x 10 12 G. From this figure, one 
can see that the effective potential is positive, therefore due to the 
positivity of the effective potential (V e // > 0) neutrinos can os¬ 
cillate resonantly. Taking into account the parameters of the two 
(table 2) and three-neutrino (table 3) mixing we analyze the reso¬ 
nance condition, as shown in Figure 5. From this figure we can see 
that temperature is a decreasing function of chemical potential, for 
the values of temperature in the range 0.7 MeV ^ T ^ 5 MeV. We 
found that the chemical potential is in the range of 1.7 eV ^ /i ^ 
32 eV for solar (left-hand panel above), 2.8 eV ^ /i ^ 91 eV for 
atmospheric (right-hand panel above), 36 keV ^ ^ 465 keV for 

accelerator (left-hand panel below) and 187 eV < fi ^ 6.5 x 10 3 
eV, for three-neutrino (right-hand panel below) mixing. It can also 
be seen from fig. 5 that the chemical potential achieves the largest 
value when accelerator parameters are considered, and the smallest 
one if solar parameters are taken into account. 

Recently, Fraija et al. (2014) showed that neutrinos can oscillate 
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Table 3. Here we show the neutrino flavor ratio in each zone of the collapsing star for E u = 5, 10, 15 and 20 MeV. 


E u 

On the NS surface 

Accretion material 

Free fall zone 

Outer layers 

On Earth 

5 

1.2:0.9:0.9 

1.176:0.912: 0.912 

1.151:0.925:0.925 

1.124:0.938:0.938 

1.039:0.978:0.984 

10 

1.2:0.9:0.9 

1.160:0.920:0.920 

1.130:0.935:0.935 

1.113:0.944: 0.944 

1.035:0.980: 0.986 

15 

1.2:0.9:0.9 

1.158:0.921:0.921 

1.132:0.934:0.934 

1.129:0.936:0.936 

1.040:0.977:0.984 

20 

1.2:0.9:0.9 

1.177: 0.911:0.911 

1.147:0.926:0.926 

1.146:0.927:0.927 

1.045:0.973:0.982 


Parameters 

Sm 2 (eV 2 ) 

angle 

Solar 

5.6±i;» x icr 5 

tan 2 0 = 0.427l°;° 3 | 

Atmospheric 

2 -llo .4 X icr 3 

sin 2 20 = 

Accelerator 

0.5eV 2 

sin 2 0 = 0.0049 


Table 1. The best fit values of two-neutrino mixing (solar, atmospheric and 
accelerator neutrino experiments) (Aharmim & et al. 2011; Abe & et al. 
2011b; Araki & et al. 2005; Shirai & KamLAND Collaboration 2007; The 
KamLAND Collaboration & Mitsui 2011). 


Parameters 

Best fit 

Sml i 

7.62 x 10 -5 eV 2 

5m 2 1 

2.55 x 10 -3 eV 2 

sin 2 #i 2 

0.320 

sin 2 #23 

0.613 

sin 2 #13 

0.0246 


Table 2. The best fit values of the three-neutrino mixing (Forero et al. 2012). 

resonantly due to the density profiles of the collapsing material sur¬ 
rounding the progenitor. The neutrino effective potential associated 
to each region is V e ff = V2 GfNa p(r)Y e with Y e the number of 
electron per nucleon and p(r) given by eqs. (10), (13) and (14). We 
plot the survival and conversion probabilities for the active-active 
(zz e ,/x,T <-> ^e,p,r) neutrino oscillations in each region as shown 
in Figure 6. From top to bottom panels, we show the oscillation 
probabilities when neutrinos pass through the new thin crust on 
the neutron star surface (top), the quasi-hydrostatic envelope (sec¬ 
ond), the free fall zone (third) and the external layers (four). One 
can see that these probabilities vary with neutrino energy (left) and 
distance (right). Taking into account the oscillation probabilities in 
each region and in the vacuum (on its path to Earth), we calcu¬ 
late the flavor ratio expected on Earth for four neutrino energies 
(E u = 5 MeV, 10 MeV, 15 MeV and 20 MeV), as shown in table 
3. In this table we can see a small deviation from the standard ratio 
flavor 1:1:1. In this calculation we take into account that for neu¬ 
trino cooling processes (electron-positron annihilation, inverse beta 
decay, nucleonic bremsstrahlung and plasmons), only inverse beta 
decay is the one producing electron neutrino. It is worth noting that 
our calculations of resonant oscillations were performed for neutri¬ 
nos instead of anti-neutrinos, due to the positivity of the neutrino 
effective potential. 

Finally, we estimate the number of initial neutrino burst expected 
during the neutron star formation. Considering a temperature T^ 4 


MeV, a duration of the neutrino pulse t ~ 10 s, the average neu¬ 
trino energy < E Pe 15 MeV and a total fluence equivalent of 
Cas A $ « 1.02 x 10 12 v e cm -2 (Mohapatra & Pal 2004; Bah- 
call 1989), then the total number of neutrinos emitted from Cas A 
would be Ntot — 6<F47rd 2 « 8.48 x 10 57 . Similarly, we can 
compute the total radiated luminosity corresponding to the binding 
energy of the neutron star L u « x < E Pe 2 x 10 52 erg/s. 

Taking into account the effective volume V ~ 0.56 x 10 12 cm 3 
(Hyper-Kamiokande Working Group et al. 2014), from eq. (16) 
we get N ev ~ 2.29 x 10 6 events that could have been expected 
during the neutron star formation in a neutrino detector as Hyper- 
Kamiokande experiment. Comparing the number of neutrinos that 
would have been expected during the neutron star formation and 
the hyper accretion phase, we would get 716.75 events. 


5 CONCLUSIONS 

Although the Cas A supernova event occurred 330 years ago, its 
similarity to other core-collapse supernovae allows to study the 
early evolution of the proto-neutron star inside the supernova rem¬ 
nant. In the present work, we have done a study of the hypercritical 
phase onto the newborn neutron star, in the Cas A scenario, in order 
to analyze the neutrino cooling effect when a quasi-hydrostatic 
envelope is formed. This phase is important because it allows 
to submerge the magnetic field within the new crust formed by 
the strong accretion. In this context, after the hyperaccretion has 
stopped, the magnetic field could diffuse back to the surface and 
result in a delayed switch-on of a pulsar (Muslimov & Page 1995). 
Depending on the amount of accreted matter, the submergence 
could be deep enough that the neutron star may appear and remain 
unmagnetized for several centuries or millennia (Geppert et al. 
1999). Cas A has been classified as one of Central Compact 
Objects (CCO; Pavlov et al. (2004)) defined as X-ray sources with 
thermal-like spectra observed close to the centers of supernova 
remnants without any counterparts in radio and gamma wavebands. 
With blackbody temperatures of hundreds of eV and luminosities 
in the range 10 33 — 10 34 erg s -1 , they present no evidence of 
a pulsar magnetic field. This scenario was recently revisited by 
Vigano & Pons (2012) in the context of diffusion of the magnetic 
field post-hyperaccretion. Following these ideas, we perform 
hydrodynamic numerical simulations of the hypercritical phase 
for the Cas A scenario, focusing on the formation of a thin new 
crust built in such a phase. We found that the size of this crust is 
at the same height-scale of the neutrino-sphere where the neutrino 
cooling is operative. We show that a very necessary ingredient for 
forming an atmosphere in quasi-hydro static equilibrium is a sink 
of energy at the bottom of the envelope, which allows to deposit 
large amounts of material onto the neutron star’s surface. The 
copious amount of neutrinos generated on the neutrino-sphere, at a 
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very efficient rate, provides this support. 

As a signature of this phase we calculate that the number of events 
expected on Hyper-Kamiokande detector is 3195. This number 
of events going through the supernova remnant will oscillate, 
firstly in the thermal and magnetized plasma, secondly, due to 
different electron densities from two to four regions and finally, on 
their paths (vacuum) to Earth. For the first region, we have used 
the neutrino effective potential derived in Fraija (2014) which is 
a function of temperature (T), chemical potential (/i), neutrino 
energy (E u ) and magnetic field ( B ). We have shown that for a 
neutrino test of energy 1, 5 , 10 and 20 MeV, and parameters 
considered of temperature and chemical potential in the range of 1 
MeV ^ T ^ 5 MeV, 1 eV ^ /i ^ 4.3 x 10 2 eV, neutrinos oscillate 
resonantly, for two- and three-neutrino mixing. In regions from 
two to four, we have also calculated each effective potential and 
then analyzed their oscillations through these regions. Considering 
the values obtained in the simulations for Cas A, the neutrino 
flavor ratios expected on Earth were computed. Our analysis shows 
that deviations from 1:1:1 are obtained for neutrino energies of 5, 
10, 15 and 20 MeV (see table 3). Diverse flavor ratios of thermal 
neutrinos will give us information about the parameters involved 
in this hypercritical accretion episode. 

RCW103, Pup A, and Kes 79 (Kaspi 2010) supernova remnants 
are good candidates to test the validity of the hypercritical phase 
and the subsequent submergence of the magnetic field, in the early 
history of these supernovae. 
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Figure 1. Pseudo-color of density (A) and neutrino emissivity (B), for the Cas A parameters: the initial transient, t = 0.2 ms (up); shock evolving in the 
domain, t = 1 ms (right); shock leaves the domain and the transient begins to disappear, t = 5 ms (down); and a quasi-hydrostatic equilibrium envelope is 
formed with a new crust on the stellar surface, t = 30 ms (left). 
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Figure 2. Density and temperature profiles as a function of the stellar radius. In the left-hand panels we show the density and temperature profiles obtained 
with MHD simulations with similar parameters of Cas A (Bernal et al. (2013)), whereas in the right-hand panels the analytical approach developed in Chevalier 
(1989) is shown. Here, [1] represents the Cas A accretion rate, and [10] and [100] are accretion rates one and two orders of magnitude larger, respectively. 
Note the new crust formed in the hypercritical regime very close to the neutron star surface (~ 300 m). The quasi-hydrostatic envelope lies between the new 
crust until ~ 50 km. The free-fall region is the remaining region, for the Cas A parameters. 
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Figure 3. Left-hand panel above: Magnetic Field Profiles as a function of radius are plotted. Note the complete submergence of the magnetic field at the same 
height scale of the neutrino-sphere. Here, [1] represents the Cas A accretion rate, and [10] and [100] are accretion rates one and two orders of magnitude larger, 
respectively. Right-hand panel above: Parameter space of temperature and density, and regions where each neutrino cooling process is dominant. The white- 
dashed lines divide the region of neutrino production for Cas A. Panel below: We show the neutrino luminosity as a function of time for the hydrodynamic 
case. This plot was obtained integrating over the entire computational domain and using the tabulated numerical values performed by Itoh et al. (1996) . 
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Figure 4. Neutrino effective potential is plotted as a function of magnetic field for three values of temperatures (1, 2.5 and 5 MeV). We have used the neutrino 
energy of 5 MeV 
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Figure 5. Contour plots of temperature and chemical potential as a function of neutrino energy for which the resonance condition is satisfied. We have required 
the best-fit values of the two-neutrino mixing (solar, top left; atmospheric, top right; and accelerator, bottom left) and three-neutrino mixing (bottom right). 
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Figure 6. Oscillation probabilities of neutrinos as a function of neutrino energy (left) and distance (right) for the three-neutrino mixing parameters. We show 
the four regions (from top to bottom panels): The new crust on the neutron star surface (top), quasi-hydrostatic envelope (second), free fall zone (third) and 
external layers (four). 
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